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ABSTRACT 


An  analytical  investigation  and  a  numerical  procedure  are  presented 
for  calculating  the  fully  viscous  shock -layer  flow  near  the  sharp  leading 
edge  of  both  two-dimensional  and  axisymmetric  slender  bodies  under 
hypersonic  low -density  flow  conditions.  The  calculation  method  is  based 
on  numerical  integration  of  the  governing  continuum  fluid  mechanical 
equations  of  motion  using  an  implicit  finite -difference  technique  in  con¬ 
junction  with  wall  slip  boundary  conditions;  the  shock  structure  is  deter¬ 
mined  as  part  of  the  resulting  solution.  The  mathematical  model  used 
in  the  present  investigation  represents  a  considerable  improvement  over 
previously  published  continuum  studies  in  that  a  more  complete  and 
accurate  set  of  governing  equations  and  boundary  conditions  is  employed. 
The  results  of  the  numerical  calculations  have  been  compared  with  pre¬ 
vious  flat-plate  test  results  taken  in  the  von  Karman  Facility  Low  Density 
Hypersonic  Wind  Tunnel  (L)  by  Becker  and  Boylan;  in  general,  good  to 
excellent  agreement  between  the  present  theory  and  experiment  is 
revealed. 
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NOMENCLATURE 


-A-n,  ®n*  Cn,  Dn 
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P 

Po 
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Re 


Coefficients  in  finite -difference  equations,  Eq.  (39) 
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Stanton  Number,  q/'pJJa>  cp  (T0  -  T  J 

Coefficients  in  inversion  formulas,  Eq.  (44) 

Static  enthalpy 
Thermal  conductivity 

Characteristic  length,  1  in.  in  all  calculations 
Molecular  mean  free  path 
Mach  number,  Um  /VtR  T. 

Total  number  of  grid  points 
Prandtl  number,  cpWk 

Static  pressure 
Pitot  pressure 
Heat  flux 
Gas  constant 

Reynolds  number  p^UJL/jUo, 


Reynolds  number  p^U^x/p^ 

Radius  in  cylindrical  coordinates 
Local  radius  of  cone 
Static  temperature 
Stagnation  temperature 
Free -stream  velocity 
x- component  of  velocity 
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Vx  a  Hypersonic  interaction  parameter,  Ma  n/C/  vRe^ 

v  y-component  of  velocity 

w  General  dependent  variable 
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tt3»  °4 
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Coordinate  along  the  body 

Coordinate  normal  to  the  body 

Thermal  accommodation  coefficient 

Coefficients  in  "standard"  parabolic  differential 
equation 

Ratio  of  specific  heats 

Transformed  normal  coordinate,  Eq.  (18) 

Angle  of  attack 
Cone  half  angle 
Dynamic  viscosity 
Bulk  viscosity 

Transformed  independent  variable 

Pi,  3.  14159 

Density 

Accommodation  coefficient  for  momentum 
Hypersonic  parameter,  \Zc/-s/Rex 
Denotes  dimensional  quantities 
Denotes  free -stream  conditions 
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SECTION  I 
INTRODUCTION 


1.1  BACKGROUND 

Rarefied  hypersonic  flow  over  slender  bodies  having  sharp  leading 
edges  has  been  a  problem  of  considerable  interest  for  a  number  of 
years.  A  great  deal  of  both  experimental  and  theoretical  work  has  been 
done  on  flows  of  this  type;  in  particular,  the  hypersonic  flow  over  a  flat 
plate  has  received  much  attention.  From  these  studies  a  fairly  definite 
qualitative  description  of  the  flow  field  has  emerged.  A  schematic 
showing  the  essential  features  of  hypersonic  flow  over  a  flat  plate  is 
presented  in  Fig.  1.  It  is  well  known  that  the  viscous  retardation  of  the 
flow  near  the  plate  surface  deflects  the  outer  streamlines  and  generates 
a  shock  wave.  The  so-called  merged  layer  begins  where  the  molecular 
mean  free  path  is  small  enough  compared  to  the  thickness  of  the  disturbed 
layer  that  the  equations  of  continuum  fluid  mechanics  are  applicable. 
Upstream  of  the  merged  layer,  near  the  leading  edge,  a  noncontinuum 
approach  is  required.  In  the  merged  layer,  the  shocklike  structure  is 
still  diffuse  and  cannot  be  treated  separately  from  the  viscous  region. 

At  the  surface  of  the  plate  there  is  relative  slip  motion  between  the  plate 
and  the  fluid  and  a  temperature  jump  across  this  interface.  As  the  flow 
proceeds  downstream,  the  compression  region  merges  into  a  Rankine- 
Hugoniot  shock,  and  a  distinct  inviscid  region  develops  between  the  shock 
wave  and  the  viscous  boundary  layer.  This  region  is  known  as  the  strong 
interaction  regime,  where  the  classical  boundary -layer  equations  are 
applicable,  but  the  development  of  the  viscous  flow  region  and  the  invis¬ 
cid  flow  are  still  strongly  interdependent.  Further  downstream  the  self- 
induced  shock  wave  weakens  and  the  inviscid  flow  is  essentially  a  per¬ 
turbation  of  the  undisturbed  flow  conditions.  In  this  weak  interaction 
region  classical  Prandtl  boundary-layer  theory,  coupled  to  the  inviscid 
flow  by  the  displacement  effect,  can  be  applied. 

The  present  study  is  concerned  with  hypersonic  flow  in  the  merged 
layer  regime,  where  a  continuum  analysis  is  applicable  but  the  Prandtl 
boundary -layer  equations  are  inadequate  to  properly  describe  the  flow. 

In  this  report  a  calculation  procedure  for  two-dimensional  hypersonic 
flow  over  flat  plates  and  wedges  and  axisymmetric  flow  past  circular 
cones  is  presented. 
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Shock  Wave 


1.2  PREVIOUS  ANALYTICAL  STUDIES 

One  of  the  first  efforts  to  analyze  the  merged  layer  over  a  flat 
plate  was  carried  out  by  Pan  and  Probstein  (Ref.  1).  The  mathe¬ 
matical  model  employed  the  usual  compressible  bound  ary -layer  equa¬ 
tions  augmented  by  the  inviscid  y-momentum  equation  to  describe  the 
transverse  pressure  variations.  These  equations  were  reduced  to 
ordinary  differential  equations  by  a  local  similarity  approximation. 
Velocity  slip  and  temperature  jump  boundary  conditions  were  applied 
at  the  plate  surface.  Taking  into  account  shock  thickness  and  curvature, 
modified  Rankine-Hugoniot  conditions  were  derived  for  the  tangential 
component  of  velocity  and  temperature  behind  the  shock.  These  condi¬ 
tions  were  used  as  the  outer  boundary  conditions  for  the  merged  layer. 

Shorenstein  and  Probstein  (Ref.  2),  following  the  work  of  Oguchi 
(Ref.  3),  extended  the  analysis  of  Pan  and  Probstein  to  take  into 
account  the  departure  of  normal  velocity  component  and  density  from  the 
Rankine-Hugoniot  values.  Shorenstein  and  Probstein  used  the  same  set 
of  simplified  equations  as  Pan  and  Probstein  to  describe  the  flow  in  the 
merged  layer.  The  results  of  both  of  these  studies  showed  good  quali¬ 
tative  agreement  between  the  data  compared.  The  quantitative  agree¬ 
ment  generally  improved  as  the  effects  of  shock  curvature  and  wall  slip 
were  included  in  the  calculations. 

Chow  (Ref.  4)  divided  the  merged  layer  into  two  subregions,  a 
shock  wavelike  region  and  a  boundary-layerlike  region,  requiring  that 
the  velocity  and  the  velocity  gradient  be  continuous  along  their  common 
boundary.  An  integral  analysis  was  applied  to  the  merged  layer,  assum¬ 
ing  a  cubic  velocity  profile  with  slip  in  the  boundary -layer  region.  The 
Crocco  relationship  was  used  to  express  temperature  in  terms  of 
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velocity  and  thereby  eliminate  the  energy  equation.  Comparison  of  the 
results  of  the  calculations  with  experimental  data  showed  reasonably 
good  agreement. 

In  a  series  of  reports  and  publications  (Refs.  5,  6,  and  7)  Rubin 
and  his  colleages  have  documented  a  calculation  method  that  is  very 
similar  to  that  of  the  present  study  and  indeed  have  provided  the  incentive 
for  this  work.  In  the  initial  publication  in  the  open  literature,  Rudman 
and  Rubin  (Ref.  6)  reduced  the  governing  conservation  equations  to  a 
form  applicable  to  the  entire  disturbed  flow  region  for  a  two-dimensional 
flat  plate.  The  equations  used  were  essentially  two-dimensional 
boundary -layer  equations  ignoring  the  longitudinal  pressure  gradient 
3p/3x  but  augmented  by  the  lateral  momentum  equation  to  account  for  the 
transverse  pressure  gradient.  These  equations  were  integrated  step  by 
step,  using  an  explicit  finite -difference  technique  through  the  entire  dis¬ 
turbed  layer  including  the  shocklike  structure.  The  results  of  these 
calculations  agree  fairly  well  with  experimental  pressure  and  heat- 
transfer  data  although  the  profiles  exhibited  an  unusual  behavior  in  the 
shock  wave  region  as  the  strong  interaction  region  was  approached. 

This  was  initially  thought  by  Rudman  and  Rubin  to  be  a  weak  compres¬ 
sion  wave  behind  the  shock  but  was  later  found  to  be  a  result  of  an  in¬ 
adequate  description  of  the  shock  structure.  In  a  succeeding  publication, 
the  theory  was  extended  to  axisymmetric  cones  and  simple  three- 
dimensional  flows  (Ref.  7). 

Cheng  et  al.  (Ref.  8)  developed  a  calculation  procedure  very  simi¬ 
lar  to  that  of  Rudman  and  Rubin  for  application  to  the  merged  layer  for  a 
flat  plate.  The  equation  set  used  was  somewhat  different  from  that  of 
Ref.  5,  however;  in  particular,  the  longitudinal  pressure  gradient  term 
was  included  in  the  x-momentum  equation.  Inclusion  of  this  term  can 
produce  a  singularity  in  the  equations  at  the  sonic  point  u  =  a.  This 
singularity  is  discussed  in  some  detail  by  Cheng  et  al.  and  by  Token 
(Ref.  9),  who  used  a  similar  equation  set  to  calculate  viscous  hyper¬ 
sonic  flow  in  slender  axisymmetric  nozzles.  The  profiles  obtained  by 
Cheng  et  al. ,  like  those  of  Rudman  and  Rubin,  exhibited  a  kink  in  the 
outer  portion  as  the  shock  structure  became  thin.  It  was  reported  that 
this  kink  disappeared  as  the  grid  and  step  sizes  of  the  calculation  pro¬ 
cedure  were  reduced. 

Recently  Kot  and  Turcotte  (Ref.  10)  presented  an  interesting  study 
on  flow  over  a  flat  plate  in  which  the  hypersonic  free  stream  is  modeled 
by  a  zero -temperature  molecular  beam.  By  collision  with  continuum 
molecules,  these  beam  molecules  are  themselves  converted  to  continuum 
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molecules.  The  beam  molecules  are  assumed  to  have  a  delta-function 
velocity  distribution  and  the  continuum  molecules  a  near-Maxwellian 
distribution.  Kot  and  Turcotte  developed  partial  differential  equations 
describing  the  interaction  of  the  beam  molecules  and  the  continuum 
fluid  and  integrated  these  equations  using  a  finite -difference  scheme. 
The  results  of  their  calculations,  however,  appear  to  be  very  similar 
to  pure  continuum  approaches  such  as  those  used  in  this  study. 


SECTION  II 

MATHEMATICAL  MODEL 


2.1  GOVERNING  EQUATIONS 

The  mathematical  formulation  will  first  be  presented  for  two- 
dimensional  hypersonic  flow  over  flat  plates.  The  geometry  and 
coordinate  system  under  consideration  for  this  case  are  shown  in 
Fig.  2. 


b.  Nonzero  Angle  of  Attack 

Fig.  2  Geometry  and  Coordinates  for  Two-Dimensional  Flow 
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The  equations  describing  such  a  flow  field  are  the  steady,  two-dimensional 
compressible  conservation  equations  for  mass,  momentum,  and  energy  as 
given  in  Ref.  11.  Since  a  forward -marching  method  of  solution,  such  as 
is  used  for  laminar  boundary  layers,  was  planned,  the  equations  were 
parabolicized  by  deleting  second-derivative  terms  in  the  streamwise  (x) 
direction.  From  previous  analyses  (Refs.  5  and  8),  it  is  well  known  that 


for  thin,  viscous  shock  layers,  terms  of  the  form 


_0 

8x 


9  /.-  8w\ 

8X7 


are  negligible. 


compared  to  other  terms  of  the  equations.  It  should  be  pointed  out  that 
under  the  thin  shock-layer  assumption  there  are  in  the  conservation  equa¬ 
tions  additional  terms  which  are  of  the  same  order  as  these  terms  and 
which  could  be  justifiably  dropped  on  an  order- of- magnitude  basis  (Ref.  5). 
However,  since  the  additional  terms  could  be  handled  by  the  technique 
used  to  solve  the  equations,  it  was  considered  worth  the  slight  extra  pro¬ 
gramming  effort  and  machine  time  required  to  retain  them.  With  the  terms 
_8_ 

8x 


9  /  ffw  \  9  /  -  9  T  \ 

—  fju  and  —  ^k  deleted,  the  governing  equations  are  as  follows: 


Continuity  Equation: 


2e1  +  =  o 

dx  dy  (1) 


x-Momentum  Equation: 


(2) 


y-Momentum  Equation: 
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Energy  Equation: 


The  fluid  under  consideration  is  assumed  to  be  a  perfect  gas  with 
the  following  constant  specific  heats: 


P  =  pET 


(5) 


dh  =  cp  dT 


(6) 


The  dynamic  viscosity,  M,  and  the  bulk  viscosity,  are  both 
assumed  to  be  functions  only  of  the  static  temperature,  T.  The  mathe¬ 
matical  model  is  completed  by  specifying  boundary  and  initial  conditions 
for  the  dependent  variables.  At  the  wall,  standard  slip  and  temperature 
jump  conditions  (Ref.  12)  applicable  to  low-density  hypersonic  flows  are 
employed  along  with  a  condition  of  zero  normal  velocity  at  the  plate  sur¬ 
face. 


u(x,0) 


T(x,0) 


2  -  a  2y  (l  <3t\ 

a  y  +  l\Pr  dy/ _ 

y=0 


+  T’.M 


v{x,0)  =  0 


(7) 

(8) 

(9) 


In  the  above  equations,  a  and  a  are  the  accommodations  coefficients  for 
momentum  and  energy,  respectively,  and  S.  is  the  molecular  mean  free 
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path.  The  accommodation  coefficients  are  assumed  to  be  unity  corre¬ 
sponding  to  diffusely  reflected  molecules  with  perfect  thermal  accommo¬ 
dation.  The  thermal  creep  term  of  Eq.  (7)  is  frequently  neglected  in 
low -density  hypersonic  flow  studies. 

Outside  the  shock  layer,  the  flow  variables  assume  their  free- 
stream  values 

*  v-0,  T-T^  ,  p-*px  (10) 

For  the  case  of  a  plate  at  angle  of  attack,  the  velocity  components  out¬ 
side  the  shock  layer  become 

u  -»  U  cos  6  and  "v  — U  sin  0 

oo  oo 

where  6  is  the  angle  between  the  plate  surface  and  the  free  stream  as 
shown  in  Fig.  2. 

The  equations  are  now  put  in  dimensionless  form  by  the  introduction 
of  the  following  dimensionless  variables: 


The  dimensionless  equations  and  boundary  conditions  that  result  are  as 
follows: 

Continuity  Equation: 


dp  u  dp  v 
dx  +  dy 


0 


x-Momentum  Equation: 


d  u  da 

P»T*  +  P'Ty  = 


1 


yK 


+ 


i 

Re 


d2u  ft  d2v  dp  d T  da 

P  fiy2  +  3  dxdy  +  dT  dy  dy 


dft  dT  dv  2  dft  dT  dv  _d_  /  dv\ 
dT  dy  dx  3  dT  dx  dy  dx  \  B  dy) 


U2) 


(13) 
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y-Momentum  Equation: 


du  du  1  /  dT  ~  dp\  1  [~4  d2v  u.  d2u  4  dii  dT  By 

pu  3x  +  pv  By  ~  yM2  V  By  +  By)  +  +  3  Bxdy  +  3  dT  By  By 

_  1  (14) 

2  d ji  dT  d_u  dp  dT  du  _d_  /  du\  d  /  dv\] 

3  dTdy  dx  +  dT  dx  dy  +  dy  \^Bxj  +  By  \  Bdy/| 


Eneigy  Equation: 


dT  dT 

paB7  +  ^vd7 


,  T/du  dv\  y  d2 T  dti  /dTV 

y  P  Vx  *  By)  +  RePr  P  <jy2  +  dT  \dy/ 


y(y  -  1)M2 

+  — r, — P 


[~4/duV  4  / dv\2  ( dv\2  /du\2  „  dv  du  4  du  <9 v~| 

[3\dx/  +  3  \dy/  +  \Bx)  +  \dy/  +  dx  dy  3  dx  dyj 


y(y  —  DM,,,,  /(Ju\  „  du  dv  /dA2^ 

Re  ^®\dx/  +  dx  dy  +  \dy/ 


(15) 


In  these  equations,  the  pressure  has  been  eliminated  in  terms  of  density 
and  temperature  by  using  the  equation  of  state. 

The  dimensionless  boundary  conditions  are 


u(x,0) 


2  -  a 
a 


y=0 


+ 


3  _1_  j_/dT\ 

4  Re  pT  \dx/y_0 


T(x,0) 


2  -  a  2y  J_f  dT\ 
a  y  -  1  PrVy/pQ 


+  Tw(x) 


(16) 
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where 


£ 


Z_ 

L 


>/y  rw  jj_ 

Re  * 2  pvt7 


(17) 


Outside  the  shock  layer  the  dependent  variables  u,  T,  and  p  all  go  to 
unity  and  v  vanishes  for  zero  angle  of  attack.  If  the  plate  is  inclined 
to  the  free  stream  at  angle  6,  then  outside  the  shock  layer 

p  -*  1  ,  T  -*  1 

u  ->  cos  6  ,  v  -*  —sin  6 


It  should  be  noted  that  the  explicit  Reynolds  number  dependence 
can  be  removed  from  the  equations  by  defining  new  dimensionless  in¬ 
dependent  variables 


x*  =  xRe 
v*  =  yRe 

The  equations  and  boundary  conditions  in  x*  and  y*  variables  are  the 
same  as  those  in  Eqs.  {12)  through  (17)  with  Re  set  equal  to  one.  Thus 
the  effect  of  varying  Reynolds  number  on  the  solution  is  simply  to 
change  the  length  scale. 


2.2  TRANSFORMATION  OF  EQUATIONS 

In  order  to  approximate  the  governing  partial  differential  equations 
accurately  with  finite-difference  equations,  it  is  necessary  that  the 
grid  point  spacing  be  very  small  in  the  compressive  shocklike  region 
where  abrupt  changes  occur.  It  is  desirable,  therefore,  to  transform 
to  a  new  set  of  independent  variables  which  will  expand  the  compres¬ 
sive  region  and  make  the  equations  more  amenable  to  a  finite-difference 
solution.  Since  the  density  is  at  a  maximum  in  the  shocklike  region,  the 
Howarth  transformation 


f  =  X  1)  =  frpiy 
0 


(18) 
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produces  the  desired  deformation.  Under  this  transformation  the 
governing  equations- take  the  following  forms: 

Continuity  Equation: 


dpu  drj  dp u 
d£  +  <9  x  drj 


dpv  n 
+  p  -k~  =  0 
dr) 


(19) 


x-Momentum  Equation: 


d»  9  <?u  ]  /  dT 

par» +  p  = 


(\  ^b\  d2v  2  dp  dT  <?u 

+  *\3  +  p/dxdy  +  P  ATdrj  drj 

/^B  2\  dp  dv  dT*] 

+  \p  3 )P  dT  drj  dxj 


dp  dT  dv 
+  PdT  d^  dl 


(20) 


y-Momentum  Equation: 


d v  9  dv 

pud7  +  '  v*i  “ 


dp  dv\ 
P  dj?/ 


/l  ^b\  dzu  /4  ^b\  2  dT  dv 

P\3  p/dxdy  +  \3  +  p)P  AT  drj  drj  (21) 

(Pb  2\  dp  dT  du  dp  dT  du"l 

+  \  p  3/P  dT  drj  dx  +  ^dT  dx  d^J 


Energy  Equation: 
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where 


d  d  dr]  d 

dx  d£  +  dx  dr/ 

d 2  _  92  dr)  d 2  d2T)  d 

dxdy  ^  d^drj  +  ^  dx  ^ 2  <9 xe?y  drj 


(23) 

(24) 


Equations  (19)  through  (22)  are  the  equations  that  were  programmed  for 
solution.  It  is  evident  that  the  Howarth  transformation  substantially 
complicates  the  equations  algebraically.  In  particular,  the  derivative 

&TI 

does  not  cancel  out  as  in  the  boundary -layer  case;  hence  this  term 
must  be  evaluated  in  the  course  of  the  solution. 


2.3  CONICAL  FLOW 

The  governing  equations  for  axially  symmetric  flow  about  a  sharp 
conical  body  have  also  been  formulated  for  solution  on  a  digital  com¬ 
puter.  The  geometry  and  nomenclature  for  flow  about  a  sharp  cone  are 
shown  in  Fig.  3. 


Fig.  3  Geometry  and  Coordinate  System  for  Axisymmetric  Flow 


The  continuity  equation  for  axially  symmetric  flow  in  the  notation  shown 
above  is 


l^(pur)  +  r)  = 


(25) 
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After  transforming  to  (1,7?)  variables  the  continuity  equation  is 


dp u  dn  dpu 
w  +  dx  Tf 


+ 


o 


(26) 


This  equation  differs  from  the  continuity  equation  for  two-dimensional 
flow,  Eq.  (19),  only  by  the  final  term.  In  like  manner  the  x-momentum 
equation,  (20),  the  y-momentum  equation,  (21),  and  the  energy  equa¬ 
tion,  (22),  for  two-dimensional  flow  can  be  adapted  to  the  flow  about  a 
cone  by  the  addition  of  the  following  terms  to  the  right-hand  side  of 
those  equations:  To  the  x-momentum  equation,  Eq.  (20),  add 


1  *\  ■  a"l3u  a  [>B  2\  •  *  dix  ai* 

^  S,B  s  +  ^cos0c^  +  v 

/Mb  2\  si“2(9c  2 M  I/Mb  2\  „  dpt  <9T 

"  IT*  “  "T-  ~  f  Sln  0cJ  ‘  ij  C0S  9*  dT  d7 

/Mb  2\  sin  0c  cos  ec  ^M  cos  0cl  /Mb  i\ 

-\T~  ?> — 7 - — \v  +  (t  +  V 


cos  0  «  -=— 


c  ^  dx\ 


(27) 


To  the  y-momentum  equation,  Eq.  (21),  add 

1  (/Mb  l\  .  .  <?u  /Mb  4\  a  dv  a  d* 

*7  jvr  +  3rm<k’,% "  \t*v  ™ + 

To  the  Energy  equation,  Eq.  (22),  add 

(si*0c  dJ  cos0c  ^T\  /sin  $c  cos0  \ 

RePrM\r  *  -7“  "  frf)  *  <>'-  ,)pT\~r  P  +  VJ 


(28) 


y(y-  1)M; 


(29) 


+  ^Mb  +  4^Singc  u  +  cos  9C  v^J 
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g 

In  the  above  expressions,  7^  represents  the  differential  operation, 
Eq.  (23). 


Equations  (19)  through  (22)  and  corresponding  equations  (including 
the  conical  flow  terms)  were  programmed  for  numerical  solution  on  a 
digital  computer.  The  solution  technique  used  is  discussed  in  detail  in 
the  following  section. 


SECTION  III 

SOLUTION  PROCEDURE 


3.1  FINITE-DIFFERENCE  EQUATIONS 

Following  Blottner  (Ref.  13)  and  Davis  (Ref.  14),  one  can  write 
Eqs.  (20)  through  (22),  the  momentum  equations  and  the  energy  equa¬ 
tion,  in  standard  form,  which  is  convenient  for  the  numerical  solution 
technique. 


<52v 

dr,2 


<9w 


*■  a0w  +  a. 


=  0 


(30) 


The  coefficients  a\,  012,  <*3,  and  04  are,  in  general,  functions  of  the 
dependent  variables  of  the  problem  and  are  therefore  unknown.  These 
coefficients  must  be  determined  by  an  iterative  procedure  and  will, 
for  the  moment,  be  assumed  to  be  known  functions  of  the  independent 
variables,  rj.  Put  in  this  form,  the  x-momentum  equation,  Eq.  (20), 
becomes 


<?2U  5u 

-  anl  «„2U 

drf 2  ul*7  u2 


+ 


S3 


<9u 

+  au4  fig  - 


0 


(31) 


where 

1  dp  ^  dp  (9  T  He  u  (5i)  Rev 

a,jl  ~  p  dr}  +  p  dT  dr,  pp  dx  p 

au2  =  0 

Re  dpi  _1_  /l  ^B\  (?2v  J_  <?T  <?v  2\^_  dv  <?T 

3,13  p2pyM2  +  p2  \3  +  p  Jdxdy  +  pp  dT  dr,  dx  +  pp  \p  3/dT  dr,  dx 

u  Re 

=  ~~PT  (32) 
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The  y-momentum  equation  becomes 

d  *  v  (9  v  d  \ 

*  “v,  ^  H  “v2v  *  av3  +  av-l^  = 

where 

1  <3p  |  «|m  dT  3  lie  /  «  <3  n\ 

avl  =  p  3*  +  j  ht  ^  "  Wp  K  +  ~P 

«v  2  =  0 

3  R  r  /  2  <3T  <3p\  1  •  +  %b/M)  a  2„ 

“''3  "  "pV'iiV  3’  +P  a’'  *  p!  ^  5'5v 

_L  3(^n/V)  ~  2  «ip  <rr  du  J_  J_  4*.  ii 

+  pp  j8  i)T  <\r)  <3x  jS^  pp  (IT  d\  dr) 

u  3  He 
flv4  =  ~  PP 

=  4  +  3(p»'p) 

The  energy  equation  becomes 


(33) 


(34) 


<32t  d  r  ... 

v  *‘"'•’57  *“T! 


dr 

ai3  '  aT»gj 


o 


(35) 


where 


]  <3p  1  ilp  (3 T  v  lie  Pr  u  K e  Pr  dl\ 

+  ~  ""  ’  ppy 


11  p  dr]  p  ilT  d»7  yp 

(y  —  IHIePr/du  <3  A 
1 1 2  ~  /in v  \(3n  *  P  diy/ 


PP>' 

(y  -  I )  P  r 


T3 


^(^)2 .  ^p-(pY .  (i^)2 ,  p^y  ^  2P^ 

3  \<3x/  3  f  \(3j?/  \<3x/  '  \3rj'  K 


<3v  dii 

dx  dr] 


(36) 


6(pn/p}  -  t  ^  ^ 
+  3  ^  dx  dr] 


‘T4 


He  Pr 
PP>' 
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Equation  (30)  can  be  approximated  by  a  finite-difference  equation  at  the 
general  point  of  the  flow  field  m  +  1,  n,  using  equally  spaced  central 
differences  for  the  r\  derivatives  and  a  two -point  backward  difference 
for  the  |  derivative.  This  is  illustrated  schematically  in  Fig.  4. 


* 

>  - .  < 

>— . < 

>  < 

* 

* 

►  , 

» 

n 

►  4 

► 

!iATi 

►  i 

- 

►  « 

► 

>  — ■  i 

>  < 

1 - < 

) _ 

®  -  1  n  m  +  1 


Fig.  4  Finite-Difference  Grid 


Thus,  let 


W  «>  \V 

m+l,n  m,n 


"ST 


o(&£) 


"m+l.n+l  "m+l,n— 1 


0(At/2) 


(37) 


<J2vv\  "m+l,n+l  2"m+l,n  +  "m+l,n— 1 

Mm+i.n  =  v 


(KAtj2) 


Substituting  into  Eq.  (30),  a  finite -difference  equation  is  obtained,  as 
follows: 


'm+l,n+l  2wm+l,n  +  "m+l,n-l  m+I  .n  ** m+ 1  ■"+ 1  Wm+l,n-l 

- - -  X  ®1  2^ 


V 


+  «rl-n-m+i,  +  ° 3*+  -  «rl-  Wm+l,nA>—  - 0 
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The  coefficients  n  are  the  values  of  a  at  m  +  1,  n  and  are  assumed 

to  be  known.  In  addition  the  values  wm  n  at  the  proceeding  station  are 
assumed  known.  This  equation  for  the  three  unknowns  wm+^  n+l* 
wm+l,  nj  wm+l,  n-1  can  written  in  the  form 


where 


An  wn— 1  + 


B  w 

n  n 


+  C  w  ,  =  D 

n  n+ 1  n 


2  <  n  <  N  -  1 


(38) 


(39) 


One  such  finite -difference  equation  is  applicable  at  each  interior  point 
n  at  ?  station  m  +  1.  Information  about  the  value  of  w  at  n  =  1  and  at 
n  =  N  is  supplied  by  the  boundary  conditions.  At  y  =  0  the  boundary 
conditions  for  u,  v,  or  T  can  each  be  expressed  in  the  general  form 


w(£o) 


(40) 


Using  the  formula 


0<V> 


this  can  be  written  in  finite-difference  form  as 


where 


A jWj  +  BjWj  +  C  j  W3  —  Dj 


3k, 


=  1  + 


2A77 


B,  =  - 


2ki 

At/ 


C,  =  V1- 


1_ 

At? 


D]  “  k2 


(41) 


(42) 
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The  outer  boundary  condition  for  all  three  variables,  u,  v,  and  T,  is  of 
the  form" 


w 


N 


D 


N 


where  Djj  is  a  known  constant.  Thus  the  w  profile  at  ?  station  m  +  1  is 
given  by  the  system  of  linear  algebraic  equations 


"A1 

B1 

Cl 

0 

• 

• 

0 

W1 

' 

A2 

B2 

C2 

0 

• 

• 

0 

w2 

d2 

0 

A3 

B3 

C3 

0 

• 

0 

w3 

°3 

0 

• 

A4 

b4 

C4 

0 

0 

X 

= 

0 

• 

. 

• 

• 

. 

0 

0 

• 

• 

• 

• 

• 

0 

0 

* 

• 

• 

• 

AN-1 

BN-1 

CN-l 

WN-1 

dn-i 

0 

. 

. 

. 

. 

0 

0 

1 

WN 

dn 

* 

_ 

(43) 


Because  of  the  special  form  of  the  coefficient  matrix  there  is  a  very 
efficient  algorithm  given  by  Richtmyer  and  Morton  (Ref.  15)  that  can 
be  used  to  solve  these  equations.  The  recursion  formulas  for  w  are 
as  follows: 


where 


E  w  ,  -  F 

n  n—  1  n 


2  <  n  <  N  —  1 


Wj  =  EjW2  +  F]w3  +  Gj 


(44) 

(45) 


E, 


C2  +  A2F1 


E  = 


F_  = 


" 

+  A2Ej 

°2 

- 

G1A2 

D2 

+ 

A2E1 

- 

-c 

n 

H 

n 

+ 

A  E  . 

n  n—  1 

D 

n 

- 

A  F’  . 

n  n—  1 

B 

n 

+ 

A  E  . 

n  n— 1 

=  3,  4,  ■  •  •  N-l 


(46) 
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The  equations  for  wn  can  be  solved  using  a  double-sweep  technique; 
marching  out  from  n  =  1  calculating  the  coefficients  En  and  Fn,  and 
then  marching  back  from  n  =  N  computing  wn. 


With  this  procedure,  the  x-  and  y-momentum  equations  and  the 
energy  equation  can  be  solved  to  obtain  u,  v,  and  T  profiles.  These 
values  can  then  be  used  in  computing  improved  coefficients  aj,  a3* 
and  0-4.  The  density,  however,  must  be  computed  from  the  continuity 
equation  using  a  different  procedure.  The  method  described  below  is 
one  of  many  possible  methods  that  could  be  employed.  Expanding  the 

continuity  equation  and  solving  for  -|^  gives 


dp  p  du  p  dr]  dn  dr\  dp  p2  d\  p\  dp 

d(  u  <5?  u  d  \  drj  d\  djj  u  dt]  u  di / 


(47) 


Assuming  that  trial  u,  v,  and  p  profiles  at  m  +  1  are  known,  the  right 
side  of  this  equation  can  be  evaluated  at  each  point  n  at  station  m  +  1. 
Then  an  improved  value  of  p  is  obtained  from  the  equation 


Note  that  the  calculation  breaks  down  if  there  is  no  slip  between  the 
wall  and  the  fluid.  Slip  boundary  conditions  were  used  in  all  calcula¬ 
tions.  From  the  density  profile  it  is  possible  to  compute  values  of 

d  77 

~  for  use  in  the  next  iteration  and  at  the  same  time  perform  the  in¬ 
verse  transformation 


\ 


/ 


x  =  i 


(49) 


The  whole  calculation  procedure  can  be  repeated  to  obtain  improved 
values  of  the  flow  variables  u,  v,  T,  and  p.  This  iteration  is  continued 
until  the  values  for  successive  calculations  are  within  prescribed  limits. 
All  other  pertinent  flow  quantities  can  be  obtained  from  these  profiles. 
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3.2  COMPUTER  PROGRAM 

A  flow  chart  showing  the  steps  of  the  calculation  procedure  and  the 
basic  subroutines  of  the  computer  program  is  shown  in  Fig.  5.  (Only 
the  essential  features  of  the  computer  program  are  shown  in  the  dia¬ 
gram).  In  this  flow  chart  the  symbol  W(J)  is  used  to  represent  all  four 
flow  variables,  u,  v,  T,  and  p.  The  symbol  WO(J)  refers  to  values  at 
the  previous  station,  and  WL(J)  refers  to  values  obtained  from  the 
last  iteration.  The  function  of  the  major  subroutines  called  in  the  pro¬ 
gram  is  described  below: 

SUBROUTINE  READIN  -  Reads  the  parameters  for  the  calculations 
from  data  cards  and  defines  the  initial  profiles  for  u,  v,  T,  and  p. 

SUBROUTINE  WALL  -  Calculates  coefficients  Ej_,  Fj,  and  in 
Eq.  (45)  from  boundary  conditions. 

SUBROUTINE  COEF  -  Calculates  a  coefficients  based  on  values  from 
the  last  iteration.  Calculates  the  finite -difference  coefficients  An,  Bn, 
Cn,  and  Dn. 

SUBROUTINE  SOLVE  -  Solves  the  finite -difference  equations,  Eqs.  (43), 
by  recursion  formulas,  Eqs.  (44)  and  (45). 

SUBROUTINE  RHO  -  Calculates  density  profile  using  Eq.  (48).  Iterates 
to  convergence. 

SUBROUTINE  ETA  -  Calculates  y,  the  distance  from  the  wall,  and-^-. 

J  3x 

SUBROUTINE  RITE  -  Computes  the  pertinent  parameters  such  as  heat- 
transfer  rate  and  skin  friction  and  prints  results. 

The  calcxflations  were  generally  begun  using  uniform  initial  profiles 
for  the  gas  dynamic  variables  described  by  25  points  across  the  shock 
layer  with  a  Arj  of  0.  02.  Additional  points  were  added  as  the  calculation 
proceeded  downstream  and  the  shock -layer  thickness  increased.  The 
maximum  number  of  points  that  could  be  used  was  501,  and  this  was 
usually  sufficient  to  describe  the  shock  layer  in  the  region  of  interest 
without  increasing  A rj.  The  small  Arj  was  required  to  obtain  sufficient 
resolution  in  the  shock  wave  region  and  resulted  in  long  computation 
times.  The  time  required  to  run  a  typical  case  was  between  one  and 
two  hours  on  the  CDC  1604 -B  computer.  The  longitudinal  step  size  A? 
was  usually  set  equal  to  Arj  initially,  but  was  multiplied  by  a  growth 
factor  of  less  than  one  percent  after  each  step. 
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Fig.  5  Flow  Chart  for  Computer  Program 
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An  iteration  procedure  was  required  to  determine  the  coefficients 
of  the  linearized  governing  equations.  The  convergence  of  this  iteration 
was  determined  by  testing  the  change  in  the  solution  between  successive 
iterations.  A  convergence  criterion  of  0.  1  percent  was  found  to  yield 
acceptable  accuracy  without  an  excessive  number  of  iterations. 

Failure  of  the  calculation  procedure  occurred  for  several  cases  at 
large  values  of  x.  The  failure  seemed  to  be  associated  with  a  decrease 
in  the  slip  velocity  at  the  wall  and  was  preceded  by  oscillations  in  8v/8r? 
near  p  =  0.  These  oscillations  increased  in  magnitude  and  finally  de¬ 
stroyed  the  calculation.  As  noted  earlier,  when  the  slip  velocity  be¬ 
comes  small,  the  wall  density  as  computed  by  Eq.  (47)  becomes  very 
questionable.  It  is  likely  that  this  is  at  least  partly  responsible  for 
the  calculation  difficulties  at  large  values  of  x.  When  calculation  fail¬ 
ure  occurred,  the  hypersonic  interaction  parameter,  X,  was  always  less 
than  0.  25,  indicating  that  the  flow  was  well  beyond  the  merged  regime. 


SECTION  IV 

CALCULATION  RESULTS 


4.1  TYPICAL  PROFILES 

Some  typical  results  of  the  calculation  procedure  described  in  the 
previous  section  are  presented  in  Figs.  6  through  13.  The  free-stream 
flow  conditions  used  for  this  calculation  correspond  to  flat-plate  data 
collected  by  Becker  and  Boylan  (Ref.  16)  in  the  von  Karman  Facility 
(VKF)  Low  Density  Hypersonic  Wind  Tunnel  (L)  and  are  tabulated  below: 

=  9.19 
p«,  =  55/*He 

Jx  =  96°K 
Re/L  =  1660/in. 

y  =  L4 

T  =  300° K 

w 

Pf  =  0.71 

The  values  for  viscosity  used  in  the  calculation  were  obtained  from 
the  Sutherland  formula  with  the  constants  appropriate  for  nitrogen. 
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Figures  6  through  9  show  the  flow  variable  profiles  and  illustrate 
quite  clearly  the  transition  from  a  merged  flow  region  where  the  shock 
is  buried  in  the  viscous  region  to  the  strong  interaction  regime  where 
a  rather  well-defined  shock  exists.  The  velocity  profiles  in  Fig.  6, 
for  example,  exhibit  a  relatively  large  slip  velocity  near  the  leading 
edge  but  otherwise  appear  well  behaved.  Further  downstream  a  point 
of  inflection  develops  that  is  associated  with  the  formation  of  the  com¬ 
pressive  shocklike  region.  This  point  in  the  velocity  profile  becomes 
more  pronounced  as  the  flow  proceeds  downstream.  The  later  profiles 
show  an  inviscid  region  beginning  to  form  between  the  shock  and  the 


» 
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0.7982 


0.6066 


0.4063 


0.2133 


0.1080 


Fig.  6  Typical  Velocity  Profiles 
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viscous  boundary  layer.  The  other  flow  variables  show  a  correspond¬ 
ing  evolution  from  the  merged  region  to  the  strong  interaction  region. 
The  pressure  profiles  are  particularly  revealing  with  regard  to  the  de¬ 
velopment  of  the  shock.  Near  the  leading  edge  there  is  a  pronounced 
pressure  variation  in  the  y -direction  throughout  the  viscous  shock 
layer.  Farther  downstream  the  pressure  variation  is  largely  confined 
to  the  outer  part  of  the  shock  layer,  and  there  is  a  substantial  portion 
of  the  flow  field  where  the  boundary -layer  approximation  8p/8y  =  0  is 
valid. 


Fig.  7  Typical  Temperature  Profiles 
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The  center  of  the  shock  structure  can  be  fairly  well  defined  by  the 
peak  in  the  density  profiles  (see  Fig.  9).  The  locus  of  these  density 
peaks  is  plotted  in  Fig.  11  to  show  the  location  of  the  shock  wave  rela¬ 
tive  to  the  plate.  The  shock  is  slightly  curved  near  the  leading  edge, 
becoming  straight  and  oriented  at  13.4  deg  to  the  plate  surface  further 
downstream.  The  temperature  and  pressure  behind  a  Rankine-Hugoniot 
oblique  shock  oriented  at  the  angle  with  Ma  =  9.  19  are  shown  for  refer¬ 
ence  in  Figs.  7  and  8. 


p 


Fig.  8  Typical  Pressure  Profiles 
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The  variations  of  gas  velocity,  temperature,  and  pressure  at  the 
plate  surface  are  shown  in  Fig.  12.  The  solutions  were  started  with 
uniform  profiles  at  x  =  0  (i.  e. ,  perfect  slip);  as  expected,  the  slip 
velocity  decreases  as  distance  along  the  plate  increases.  The  pressure 
increases  rapidly  from  the  initial  free -stream  value  to  a  maximum  of 
about  6  Pa,  and  then  slowly  decreases.  The  surface  gas  temperature 
exhibits  a  similar  behavior.  Figure  13  shows  that  the  skin  friction 
coefficient  and  the  Stanton  number  display  the  expected  trends. 


Fig.  9  Typical  Density  Profiles 


25 


Fig.  11  Locus  of  Density  Peaks 
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Fig.  13  Variation  of  Skin-Friction  and  Heat-Transfer  Coefficients 


27 


AE  DC-T  R-73-1 02 


Figures  14  through  17  show  velocity,  temperature,  pressure,  and 
density  profiles  on  a  flat  plate,  a  10-deg  wedge,  and  a  10-deg  cone  under 
identical  free-stream  conditions.  The  gas  in  these  calculations  is  again 
nitrogen,  and  the  free-stream  conditions  are  tabulated  below: 


\1  =  10.15 

©o 

p»  *  2(VUg 
T=c  =  144°K 
Re/L  =  SSS^in. 


y  =  1.1 
l\v  =  300°K 

Pr  =  0.71 


Fig.  14  Comparison  of  Velocity  Profiles  for 
Flat  Plate,  Wedge,  and  Cone 
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It  is  of  interest  first  to  compare  the  flat-plate  profiles  for  this  case 
with  the  higher -density  conditions  =  9.  19  and  Re/L  =  1660/in.  shown 
in  Figs.  6  through  9.  Comparing  the  density  profiles  at  x  «  0.  8  shows 
that  the  center  of  the  shock  as  indicated  by  the  density  peak  is  approxi¬ 
mately  the  same  for  both  cases  but  that  the  shock  structure  is  not 
nearly  so  sharply  defined  for  the  lower-density  case.  For  this  low 
Reynolds  number  flow  the  shock  is  still  very  diffuse  and  buried  within 
the  viscous  boundary -layer  region. 


Fig.  14  Concluded 
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Fig.  16  Comparison  of  Pressure  Profiles  for  Flat  Plate,  Wedge,  and  Cone 
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Fig.  17  Comparison  of  Density  Profiles  for  Flat  Plate,  Wedge,  and  Cone 
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These  figures  also  show  that  shock  structure  on  the  10 -deg  wedge 
is  much  stronger  and  develops  much  more  rapidly  than  on  the  flat  plate. 
The  pressure  level  is  also  much  higher  on  the  wedge,  resulting  in  a 
higher  density  and  a  lower  slip  velocity.  While  the  thickness  of  the 
shock  layer  is  approximately  the  same  for  the  wedge  flow  and  the  flat 
plate,  the  thickness  for  the  cone  is  considerably  reduced.  This  is  typi¬ 
cal  of  conical  flows  and  is  a  result  of  the  three-dimensional  character 
of  the  flow  field. 


4.2  COMPARISON  WITH  DATA 

A  comparison  of  the  results  of  the  theoretical  calculation  method 
presented  in  this  report  with  experimental  data  from  several  sources 
was  carried  out.  These  comparisons  are  shown  in  Figs.  18  through  27. 

Figures  18  and  19  are  comparisons  of  computed  surface  pressure 
with  pressure  data  collected  by  Becker  and  Boylan  (Ref.  16)  in  Tunnel  L. 
Two  calculated  curves  are  shown  in  each  of  these  figures,  one  labeled 
"pressure  slip"  and  the  other  labeled  "no  pressure  slip.  "  The  curve 
labeled  "no  pressure  slip"  gives  simply  the  gas  pressure  evaluated  at 
y  =  0.  The  "pressure  slip"  curves  on  these  two  figures  have  been  ad¬ 
justed  according  to  the  correction  given  by  Patterson  (Ref.  17),  as 
follows: 

v=o  +  16  \T  *y/y=o  (50) 

where 

gas  pressure  at  y  =  0 
wall  pressure 
8  _ 

-  RT.  average  molecular  speed 

7 T 

molecular  mean  free  path 


p(0)  = 


F  = 


p(0)  _  .  5  T  /<?u  _  9  dA 

p  6  c  V?x  ”  dy) 


The  above  correction  was  obtained  by  assuming  the  diffuse  reflection  of 
gas  molecules  impinging  on  the  wall.  The  reflected  molecules  are 
assumed  to  be  in  equilibrium  at  the  temperature  of  the  wall.  By  con¬ 
sidering  the  transfer  of  tangential  momentum  and  translational  energy 
at  the  gas/solid  interface,  Patterson  derived  the  velocity  slip  and  tem¬ 
perature  jump  boundary  conditions  corresponding  to  Eqs.  (7)  and  (8) 
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with  the  coefficients  a  and  a  both  equal  to  1.  A  similar  treatment  of  the 
normal  component  of  momentum  across  the  interface  yields  the  relation 
between  p(0)  and  pw  given  in  Eq.  (50).  This  correction  is  therefore  con¬ 
sistent  with  the  slip  boundary  conditions  applied  to  this  study.  The  pres¬ 
sure  correction  is  applied  after  the  solution  is  obtained  and  does  not 
affect  the  calculation  of  quantities  in  the  flow  field. 

Written  in  dimensionless  form  in  terms  of  (|,rj)  variables,  Eq.  (50) 
becomes 


As  can  be  seen  from  Figs.  18  and  19,  the  pressure  slip  correction  is  a 
sizable  fraction  of  the  calculated  surface  pressure  values.  These  fig¬ 
ures  also  show  that  the  corrected  values  agree  much  better  with  the 
pressure  data  obtained  by  Becker  and  Boylan  than  do  the  uncorrected 
values.  To  the  author's  knowledge  this  is  the  first  time  that  the  signifi¬ 
cance  of  the  pressure  slip  correction  for  low -density  hypersonic  flow 
has  been  demonstrated. 


Fig.  18  Comparison  of  Calculated  Surface  Pressure  with 
Experimental  Data:  High  Reynolds  Number 
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Fig.  19  Comparison  of  Calculated  Surface  Pressure  with 
Experimental  Data:  Low  Reynolds  Number 


Figures  20  and  21  show  a  comparison  of  calculated  heat -transfer 
rate  with  experimental  data  obtained  in  Tunnel  L  by  Boylan  (Ref.  18). 
The  heat-transfer  rate  was  computed  from  the  formula 


.  /r  <?T - <?u\ 

q  -  -  V  57 


including  the  sliding  friction  term  given  by  Maslen  (Ref.  19).  Note  that 
Fig.  20  is  for  a  flat  plate  at  zero  angle  of  attack  and  Fig.  21  is  for  a 
plate  inclined  at  -10  deg  to  the  flow  stream.  In  all  cases  the  agreement 
between  calculated  values  and  the  experimental  data  is  satisfactory. 
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Figure  22  shows  the  location  of  the  center  of  shock  structure  as  de¬ 
termined  by  the  maximum  pitot  pressure  in  the  shock  layer.  The 
agreement  between  the  prediction  of  the  present  theory  and  the  experi¬ 
mental  data  of  Becker  and  Boylan  (Ref.  16)  is  seen  to  be  excellent. 


Fig.  22  Comparison  of  Calculated  Shock  Location  with 
Experimental  Data 


Figure  23  shows  a  comparison  of  computed  total  pressure  profiles 
at  three  stations  with  data  of  Ref.  16.  The  level  of  the  pitot  pressure 
peak  as  well  as  its  height  is  predicted  very  well.  These  profiles  show, 
however,  that  the  computed  thickness  of  the  shock  structure  is  con¬ 
siderably  smaller  than  that  which  actually  occurs.  The  computed  total 
pressure  profiles  for  argon  shown  in  Fig.  24  likewise  exhibit  a  thinner 
shock  structure  than  the  experimental  profiles  measured  by  Becker 
(Ref.  20).  The  location  and  level  of  the  total  pressure  peaks  do  not 
agree  as  well  with  the  data  for  this  case.  This  calculation  was  made 
for  the  monatomic  gas  to  determine  whether  rotational  nonequilibrium 
effects  could  be  responsible  for  the  thickened  shock  structure.  The 
results  of  the  comparison  would  appear  to  deny  this  explanation.  This 
underprediction  of  the  shock  structure  thickness  was  also  obtained  by 
Cheng  et  al.  (Ref.  8)  and  appears  to  be  characteristic  of  shock  models 
based  on  the  Navier-Stokes  equations. 
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Figure  25  shows  a  comparison  of  computed  density  profiles  with 
some  experimental  profiles  obtained  by  McCrosky  (Ref.  21)  in  nitrogen 
flow  at  very  high  Mach  number.  The  location  and  level  of  the  density 
peak  are  well  predicted  by  the  present  calculation  procedure.  Also 
shown  in  Fig.  25  as  solid  symbols  are  some  points  taken  from  density 
profiles  computed  by  Cheng  et  al.  (Ref.  8)  under  identical  conditions. 
Their  results  are  very  close  to  the  profiles  obtained  in  the  present 
study,  except  near  the  center  of  the  shock  structure.  In  this  region, 
Cheng  et  al.  obtained  a  density  level  somewhat  higher  than  that  calcu¬ 
lated  in  this  investigation  or  shown  in  the  data  of  McCrosky.  At  x  =  2  in. 
Cheng's  profiles  exhibit  the  spurious  kink  in  the  shock  region;  there¬ 
fore,  this  peak  is  not  plotted  in  Fig.  25.  The  better  agreement  obtained 
in  the  present  study  is  attributed  to  the  more  accurate  set  of  equations 
used  and  the  better  finite-difference  description  in  the  shock  region  due 
to  the  stretching  transformation  (Eq.  18). 


McCrosky  Data  Cheng  Calculation 
x  (Ref.  21)  (Ref.  8) 


0.3  □  ■ 

1.0  O  • 


P 


Fig.  25  Comparison  of  Calculated  Density  Profiles  with 
Experimental  Data  for  Nitrogen 
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Figures  26  and  27  show  a  correlation  of  surface  pressure  and  heat- 
transfer  data  for  flat  plates  from  several  sources,  A  replot  of  some  of 
the  data  from  the  VKF  tunnels  is  included  in  these  figures.  The  surface 
pressure  and  heat -transfer  coefficients  are  plotted  in  terms  of  the  hyper¬ 
sonic  interaction  parameter  Vx  Note  that  VXj  «  x-^^  so 

_  yj^ex 

that  large  values  of  Vx  m  occur  near  the  leading  edge  of  the  plate.  The 
hypersonic  parameter  x  is  defined  as 


X  = 


\/R 


2 

oo 


where  C  is  the  Chapman-Rubesin  constant, 

®  □  Becker  (Ref.  20),  A r 

-  11.18,  Re/E  -  244/in, 

@  O  Becker  and  Boy Ian  (Ref.  IS),  K2 

-  10.15,  Re/E  -  388/in. 

@  A  McCrosky  (Ref.  21),  »2 

8„  -  24.5,  Re/E  -  10,000/in. 

®  O  Vidal  and  Barta  (Ref.  22),  »2 

-  19.2,  Re/E  -  670/in. 


Fig.  26  Comparison  of  Calculated  Surface  Pressure  Parameter  with 
Experimental  Data  from  Several  Sources 
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^  O  Boylan  (Private  Communication),  -  18.15, 
Re/E  =  1250/in.,  Tw/Tq  -  0.1 

@  d  Same  as  above,  Small  Plate 

(5)  A  Nagamatsu  et  al .  (Ref.  23),  -  19.2, 

Re/E  =  5050/in.,  TW/TQ  =  0.235 

@  O  Vidal  and  Bartz  (Ref.  22),  =  19.2, 

Re/E  =  1350/in.,  Tw/TQ  -  0.072 

(5)  □  Vidal  and  Bartz  (Ref.  22),  *  19.2, 


Fig.  27  Comparison  of  Calculated  Heat-Transfer  Coefficient  with 
Experimental  Data  from  Several  Sources 
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The  conditions  used  to  calculate  the  solid  curves  in  these  two  fig¬ 
ures  are  designated  by  the  numbered  label  and  are  given  by  the  like- 
numbered  experimental  conditions  in  the  figure  legends. 

The  agreement  between  the  computed  surface  pressure  and  the 
experimental  data  shown  in  Fig.  26  is  considered  to  be  reasonably 
good.  The  computed  values,  however,  tend  to  be  somewhat  low  near 
the  leading  edge  of  the  plate.  The  surface  pressure  distributions  shown 
in  Figs.  18  and  19  suggest  that  this  difficulty  may  have  been  caused  by 
the  uniform  initial  profiles  that  were  used  to  start  all  calculation.  How¬ 
ever,  the  calculations  for  =  10.  15,  Re/L  =  388/in.  were  repeated 
using  different  initial  velocity  and  temperature  profiles,  and  the  surface 
pressure  distribution  was  not  significantly  altered. 

The  comparison  with  heat-transfer  data  for  flat  plates  shown  in 
Fig.  27  is  less  satisfactory.  In  particular,  the  data  of  Vidal  and  Bartz 
(Ref.  22)  is  far  below  the  corresponding  calculations.  In  view  of  the 
good  agreement  with  the  heat-transfer  datataken  inTunnelL  as  shown  in 
Figs.  20  and  21,  the  author  considers  this  to  be  an  indication  that  the 
data  reported  in  Ref.  22  may  be  somewhat  low. 


SECTION  V 

SUMMARY  AND  CONCLUSIONS 


An  analytical  investigation  has  been  carried  out  and  a  numerical 
procedure  developed  for  calculating  fully  viscous  shock -layer  flow  near 
the  leading  edge  of  two-dimensional  and  axisymmetric  slender  bodies 
under  low -density  hypersonic  conditions.  The  particular  geometries 
considered  were  two-dimensional  flat  plates  at  zero  and  negative  angles 
of  attack  and  axisymmetric  cones  at  zero  angle  of  attack.  The  calcu¬ 
lation  method  is  based  on  numerical  integration  of  the  continuum  fluid 
mechanical  equations  of  motion  using  an  implicit  finite -difference 
marching  technique.  The  mathematical  model  used  in  this  investiga¬ 
tion  represents  an  improvement  over  previously  published  continuum 
studies  in  that  a  more  complete  and  accurate  set  of  governing  equa¬ 
tions  is  employed.  Standard  slip  and  temperature  jump  boundary  con¬ 
ditions  are  employed  at  the  wall,  and  the  shock  structure  is  deter¬ 
mined  as  part  of  the  resulting  solution.  In  order  to  obtain  adequate 
resolution  in  the  shock  wave  region,  the  normal  coordinate  is  stretched 
by  using  a  Howarth  transformation. 


43 


AEDC-TR-7  3-102 


The  results  of  the  calculations  show  quite  clearly  (1)  the  develop¬ 
ment  of  a  merged  flow  regime  where  the  shock  is  buried  in  the  viscous 
region  and  (2)  the  transition  to  the  strong  interaction  regime  where  a 
well-defined  shock  structure  exists.  Near  the  leading  edge  there  is  a 
pronounced  pressure  variation  throughout  the  shock  layer,  whereas 
further  downstream  there  is  a  substantial  portion  of  the  flow  field 

r\ 

where  =  0. 

9y 

The  results  of  the  numerical  calculations  have  been  compared  with 
experimental  data  from  a  number  of  sources,  and  qualitative  agree¬ 
ment  with  surface-pressure  and  heat -transfer  data  as  well  as  with  pro¬ 
file  data  was  obtained  for  all  cases  tested.  The  quantitative  agreement 
of  the  surface-pressure  data  was  markedly  improved  by  the  application 
of  a  pressure  slip  correction  derived  by  Patterson  (Ref.  17).  A  com¬ 
parison  of  the  results  of  the  present  theory  with  experimental  data  ob¬ 
tained  on  flat  plate  and  wedges  in  Tunnel  L  reveals  very  good  agree¬ 
ment. 
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